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1. Introduction, 

The lunar Seismic Experiment Natural Activity Study has 
provided a unique opportunity to study the high frequency (4-20 Hz) 
portion of the seismic spectrum on the moon. The purpose of this 
project was to study the data obtained from the I^PE and determine 
as much as possible about l) the origin and importance of the process 
that generates thermal moonquakesj 2) the characteristics of the 
seismic scattering zone at the lunar surface; and 3) the meteoroid 
flux for masses between 0.1 and about 10 grams. 

The first problem has received the most effort and is the 
subject of the paper enclosed with this report (Appendix A). 

The detection of thermal moonquakes by the L5PE array made it possible 
to locate the sources of many events and determine that they are 
definitely not generated by astronaut activities but are the result 
of a natural process on the moon. 

The second problem, that of the propagation of seismic waves 
In the near-surface layers, was studied in a qualitative manner. 

In the absence of an adequate theoretical model for the propagation 
of seismic waves in the moon, it is not possible to assign a depth 
for the scattering layer. The I^PE data does define several parameters 
which must be satisfied by any model developed in the future. 

The third problem, that of meteoroid flux, was not studied. 

The vast numbers of events observed are almost certainly about 99^ 
thermal moonquakes and, unlike the Apollo 14 KE case (Duennebler 
and Sutton, 1974), there was no obvious way to discriminate between 
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thernal moonquakes and meteoroid impacts. At station 14, thermal 
■oonquakes always had a lower frequency spectrum than near-by 
Impacts, At the KiPE array, thermal moonquakes are found with 
all ranges of spectra. 

In addition to the above problems, a search was made for 
events detected by the PJE that were also detected oy the LSPE, 
During the times when both instruments were operatir^, only one 
event was detected that was visible on both systems. This event, 
a fairly large meteoroid Impact, was just barely visible on the 
ISPE records although easily visible at all PSE sta*.lons. The 
reason for the difference in amplitude is the different frequency 
responses of the two systems, the LJPE being less sensitive to 
frequencies below 4 Hz than the PSE, No high frequency teleselsmlc 
events (HPT) (Nakamura, et were observes during the 

times of ISPS operation, HFT events should be clearly recorded 
on the I3PE and the KE. The I5PE is now being operated solely 
for the purpose of studying events that are recorded on both the 
ISPE and PSE network, 

II, Data Collection, 

During the short periods of time curing which -.he ISPS was 
activated for the active portion of the experiment, several natural 
seismic events were noticed (Kovach, et , 1973). In the hopes 
that these events were caused by thermal moonquakes of the type 
observed at the PSE stations (Duennebler and Suttor.. 1974), the 
LSPE was activated for a 4-day period during the susvmer of 1973 
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to see If more of this activity could be observed. The results of 
this experiment showed that large numbers of thermal moonquakes 
were being detected by the I2PE. Because of the possibility of 
locating the source of many of these events using array techniques, 

It was decided to collect at least one full lunar month (lunation) 
of liiPE data. The data rates used by the liJPE require that all 
other experiments at the Apollo 1? site lose data when the KiPE 
Is activated. To prevent long periods of data loss to other exper- 
iments, the LSPE was activated for 4-day periods at different parts 
of the lunation for a period of about one year until ten 4-day 
listening i«:rlods had been obtained. In this way data covering one 
full lunation (plus two periods of overlap) were obtained. In 
addition, a 4-hour period spanning the occurrence of a solar eclipse 
was recorded. For the more than 40 days of coverage, less than 
of the data was lost due to poor transmission from the moon 
since the listening periods weie planned by JSC personnel to coincide 
with times of optimal transmission. 

Upon receipt of the data, about one full-reel tape per 4-hour 
period, compressed time scale playouts were produced (Figure l) 
using the University of Texas POP- 15 computer. These records 
were scanned by eye and events above a given size (7 mm) were recorded, 
A listing of these events and copies of the compressed playouts 
were sent to R.L. Kovach at Stanford for his use as soon as they 
became available. Once events of Interest were chosen, a set of 
event tapes was generated containing only periods where events 
occurred. These tapes were then used to generate expanded time 
scale playouts (Figures 2A and 2B of Appendix A) and power spectral 
analyses (Figure 2) of each event. 



Ill, Analysis. 
Thermal moonquakes. 
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The most time-consuming operation was the hand correlation 
of the more than one thousand events recorded, A property of 
thermal moonquakes Is that their source mechanism Is so repeatable 
that the same events are seen at about the same time of each lunation. 
Some types of events occur several times per lunation. When two 
events were found to have Identical or nearly Identical wave forms, 
they were assigned a number corresponding approximately to the 
longitude of the sun at the time of observation. Events generated 
by the LM were also identified, A complete listing of events 
observed is contained as Appendix B of this report. Of the 1730 
events listed, 324 are caused by thermal noises In the LM, and 359 
are thermal moonquakes with at least two events with nearly Identical 
character observed. It Is believed that almost all of the remaining 
events are thermal moonquakes for which no matching event was found 
either because of low slgnal-to-nolse ratio or because the period 
of recording was not long enough to record a matching event. The 
reader is referred to Appendix A for further discussion as to the 
methods used In the study of thermal moonquakes at the I2PE site. 

Seismic Propagation and Scattering , 

The eight explosive package (EP) signals and the LM impact 
signal were analyzed to obtain Information concerning the propagation 
of seismic waves and the effects of scattering. A problem encountered 
with the EP records is the occurrence of fairly high amplitude 
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electronic noise at times vfhen the detonation timing circuits are 
activated. This noise severely reduces the slgnal-to-nolse ratio. 

Cooper, 1975. used prediction filtering to remove this noise. 

Close scrutiny of the signals reveals, however, that the noise 
repeats exactly with a periodicity of 21 data points and car thus 
be removed by simple subtraction of the noise waveform from the 
signal. Power spect;~a of the corrected EP data were then calculated and 
studies of the amplitude decay with frequency and range were made. 

The results were then compared with the predictions made by Nakamura, 
et ^, , u fitf. using PSE data. His results predict that a two- 
dimensional scattering model should fit the data to ranges of at 
least 4 km. As shown by Cooper (1975) *^nd In Figure 3, this Is 
not the case. The two-dimensional scattering model Is valid to 
ranges of only one km, beyond which the observed signals reach their 
peak amplitudes well before the model would predict. The early 
peak Is apparently due to the Increase In energy from seismic waves 
passing below the scattering zone, thus a change In dlffuslvlty 
with depth that is not modelled In the two-dimensional case. 

The variation of amplitude of the signal peak with range 
also does not follow the two-dimensional scattering theory (see 
Appendix A) and is strongly frequency-dependent In such a way as 
to imply that cylindrical spreading cannot account for the decay \ 

of amplitudes beyond a range that decreases with frequency, l.e. 
high frequencies attenuate faster than expected if values of 
are used that fit the data at short range. 



6 


When this information Is coupled with that obtained from earlier 
studies, several qualitative observations concerning propagation 
of seismic waves In the moon can be made, 

1) The high degree of scattering In the surface layers prevents 
the propagation of plane waves In any direction for nore than about 

a wave length. The coherence of seismic waves on tr.e three components 
of the KE Is very low and no coherence Is observed in wave trains 
observed at the LSPE even though the separation between the closest 
sensors Is about 50 meters. These observations Imply that no 
azimuthal Information as to the source of an event can be obtained 
from observing phase variations on orthogonal Instrujients or across 
an array, 

2) The depth of the scattering zone Is very shallow (probably 
less than about 500 m and most likely less than 100 m In most regions). 
This is verified by the observation that the two-dimensional scattering 
theory does not fit the data beyond distances that decrease with 
frequency, 

3) The lunar surface Is a very poor reflector of seismic 
energy. The process of scattering coupled with the strong Increase 
In seismic velocity with depth cause almost all seismic energy tnat 
originates at or reaches the lunar surface from below to be trapped. 
Only seismic rays that leave the surface at angles within about 

one degree of vertical ever reach the lunar Interior because the 
rapid Increase In velocity with depth turns the energy back towards 
the surface. Seismic energy reaching the surface nearly vertically 
from below (such as from a deep moonquake) Is severely scattered 
and very little energy Is reflected at the near- vert leal angles 
necessary to return to the Interior, Once trapped at the surface. 
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the seismic energy decays because of absorption and slow leakage 
into the interior. This property of trapping is observed in several 
Mays, a) In general, seismic signal amplitude decreases monotonically 
with distance from the source at all times | even tens of seconds 
after an event, *he source region will contain more seismic energy 
than regions farther away, b) Events originating at the lunar 
surface have long rise times while those originating below the 
scattering layer have sharp rise times. Events (such as Impacts) 
with sources at the surface Insonlfy a region with seismic energy 
that slowly spreads across the surface of the moon and, as the 
Insonifled surface gets larger, leaks energy into tne interior. 

Events originating below the scattering layer generate silsmlc 
waves that Insonlfy a large region of the lunar surface at once 
and have short rise times, c) For events from a surface source 
high frequencies are observed to be stronger earlier in the wave 
train than low frequencies. This is because the higher frequencies 
are more easily scattered out of the surface wave guide than low 
frequencies, thus arriving at the receiver earlier. This effect 
causes a statistical dispersion effect in lunar seismograms, 
d) Surface reflections are not observed. If the lunar surface 
were a good reflector of seismic waves, then seismic phases such as 
PP, PPP, SS, and SSS would be observed. The only reliable candidates 
for these cases are at frequencies of 0.2 Hz and lower, which may 
be low enough such that the scattering layer is no longer important. 

At higher frequencies, however, the lunar surface is far too good 
an absorber to allow such phases to be seen. 
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Three factors are necessary to explain the characteristics 
of all lunar events i l) high Q in the surface regions, 2) a steep 
velocity gradient, and 3) scattering near the lunar surface. 

Without any one of these three factors the lunar seismograms would 
look vastly different. Unfortunately a model which Incorporates 
all three of these factors has not been found. Models that Include 
high Q and scattering do not Include the change in scattering with 
depth or the Increase of seismic velocity with depth. Models that 
include the Increase in velocity with depth don't Include scattering. 
Until an adequate model can be obtained it will be hard to quantify 
the above conclusions or accurately predict the characteristics 
of lunar seismograms and anelastlc properties of the r.ear-surface layers. 


Meteoroid flux . As mentioned earlier, no new information has 
been obtained on the meteoroid flux because of the lack of features 
with which to discriminate impacts from moonquakes recorded on 
the liiPE. The seismic method Itself (Duennebler, et ^. , 1975) 

Is hampered by a lack of events with which to relate amplitude 
of a seismic event generated by a meteoroid impact to the energy 
of the impacting body. The SIVB impacts, used to obtain this parameter 
In the paper above, may not accurately reflect the same energy 
partitioning as a meteoroid impact since the density of an BIVB 
is about 0,01 g/cm^ compared to a meteoroid density of greater than 
1 g/cm^. The SIVB, being less dense, will probably transfer more 
of its energy into seismic waves and less into crater formation than 
a meteoroid impact with equivalent energy. The net result would 
be an underestimate of the meteoroid flux. Since no neani. to calibrate 
this effect are available, the degree of underestimation, if any. 

Is not known. 
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IV. CONCLUSIONS. 

Th« ISPS Natural Activity Study has supplied Important Information 
concerning l) the source mechanism and Importance of thermal 
moonquakes In the degradation of the lunar surface and Z) the 
processes of scattering and propagation of seismic waves In the 
moon. It has been shown that thermal moonquakes are most likely 
caused by the movement of regollth In response to diurnal stresses 
at the lunar surface. These events would tend to move material 
from hot regions to cold regions and always In a down-slope direction, 
thus degrading lunar slopes. The rate of slope degradation Is 
still uncertain but It could be fast enough to account for a 
significant amount of the observed smoothing of the lunar surface. 

Data obtained from the explosive packages and LM Impact 

supply valuable constraints on any forthcoming models for the 
propagation of seismic waves In the moon. The need for theoretical 
work In this area must be stressed. 

There are several problems still to be solved In lunar seismology, 
and I have two main recommendations for future seismic experiments 
on the moon. 

l) Low frequency array. Because of size and weight restrictions, 
the Apollo seismometers have limited response, and seismic background 
noise at frequencies below about 1 Hz Is certainly below the threshold 
of the Instruments, Because of the very low noise conditions on 
the moon. It may be possible to detect surface waves and normal 
modes of oscillation from very small events. These Instruments 
should contain 3 orthogonal components and be spread as thickly 
and widely as possible on the lunar surface. The Information obtained 
from such an array would yield considerably more Information about 
the structure of the noon. 
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2) High frequency array. It la well known froa the Apollo 
eelsalc experiments that a large portion of the seismic energy froa 
natural lunar events Is found In the frequency band from 0,4 to 10 Hz, 
However, even in this frequency band, the Apollo Instruments could 
not detect the seismic noise level for long periods of time 
(especially during the lunar night). Thus more sensitive Instruments 
are also needed In this frequency band. 

At frequencies above 0,5 Hz there Is very little coherence 
between orthogonal components, except possibly for the first arrivals. 
In fact, little information at all Is obtained by having more than 
one component (other than redundancy), Kor this reason it would 
probably be more advisable to use the weight and data allocations 
on one high gain (vertical) Instrument at each station. To obtain 
azimuth Information it would be advisable to deploy sub-arrays 
with separations of from 2 to 10 km between components. In this way 
the phase velocity across the array could be used to obtain azimuths 
to teleselsmlc events and source locations of local events could 
be found. Efficient triggering systems could be employed to compress 
the amount of data returned, 

A system of several (5 or 6) arrays containing one set of 
orthogonal high sensitivity long period seismometers and a sub-array 
of three or four higher frequency vertical Instruments could vastly 
Increase our knowledge of the deeper Interior because of increased 
accuracy of source locations and lower threshold for event detection. 
In addition, much more can be learned about the meteoroid flux In 
the smaller mass range and about thermal moonquakes. 
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FIGURE CAPTIONS 

1 , Compressed scale playout of I*>PE data, UJPE data tapes were 
played out on a Versatec matrix plotter using the Karine 
Science Institute's PDP-I 5 computer. The numbers on the left 
correspond to the geophone numbers. The data was expanded 
from its log-compressed state and then rectified. The year, 
day, hour, and minute are labeled every 10 minutes. The 
event labeled *'KQ" is a known thermal moonquake, (The year 
should be "73" — not "72" as labeled). 

2, Spectrograms of I*SPE events. Spectrograms of three common 
types of events are shown. The top section of each figure 
shows two band-pass envelope recordings of the signal, the 
dark squares at a frequency of 19.2 Hz and the light squares 
at 9.2 Hz, The bottom section shows the power level contoured 
in 10 db steps plotted with frequency as the vertical axis 
and time as the horizontal axis. The time in days, hour, 

and minute is given at the bottom. These spectrums are not 
corrected for the frequency response of the instrument. 

Note that, compared to the thermal moonquakes, the LM event 
(center) has a sharp rise time and heavily Landed spectrum, 

3, Comparison of explosive package (EP) signal envelopes to theoreti 
cal scattering model. The theory, given by Nakamura et al. f; 
suggests that the signal amplitude should obey the function 
shown by the dashed line. While the fit is excellent at ranges 
(given in km by the number after the SP number) of less than 

1 km, the theory obviously breaks down beyond 1 km. At lower 
frequencies, the theory fits to greater ranges. 
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Thermal movement of the regolith 

Frederick Duennehikr 

The Marine Science Institute. University of Texas. Galveston. 1'exas 77^.SO 

Abstract — High frequency seismic events observed at the Apollo landing sites indicate that the repolith 
moves in response to diurnal variations of thermal stresses. While these events (thermal moonquakes) 
are observed by each of the Passive Seismic Experiment stations, the location of their sources was not 
possible until they were observed by the Lunar Surface Profiling Experiment (LSPEi array at the 
Apollo I7 site. Comparison of source locations with topographic features near the LSPE array indicate 
that thermal moonquakes are associated to some de;.,ree with large rocks and to a larger degree with 
craters. 1 hey are not associated with man-made disturbances of the lunar surface It is suggested that 
thermal motmquakes are the seismic expression of a phenomenon that is actively degrading slopes on 
the lunar surface. » 


Introduction 

TfiERMAi. MOONQU.AKES RECORDED by the ApoMo Passive Seismic Experiments 
have been studied by Duennebier and Sutton (1974). Their paper provides a 
description of the important features of the observed seismic signals including 
periodicities, spectral content, and amplitudes. They suggest that thermal moon- 
quakes are generated by fracturing of rocks along zones of weakness and by 
downslope movement of the regolith in response to thermal stresses as the surface 
temperature changes during each lunation. Because of their very small size, each 
thermal moonquake is recorded at only one station of the Passive Seismic 
Experiment, leaving no possibility of finding their source location. Estimates of 
the source distance were obtained by using knowledge of how the signal shape 
changes with distance, but this information yields only a rough guess at range, and 
no reliable correlations with surface topography were possible. 

With the deployment of the LSPE array it became feasible to locate nearby 
thermal moonquakes by standard triangulation techniques. The array consists of 
four geophones arranged in a triangle about 100 m on a side with one geophone in 
the middle. Although the prime experiment for the array was profiling of the 
surface layers using explosive sources and seismic refraction techniques (Kovach 
et ai, 1973), a four-day period of passive listening during the summer of 1973 
showed that large numbers of thermal moonquakes were being detected by the 
array. To capitalize on this fact, a series of seven addition.il four-day listening 
periods were schedulejl such that data from a complete lunation (291. dav.v) could 
be obtained over a period of about one year beginning in July. 1974. The livtoning 
periods were staggered to create a minimum of interference to other experiments 
at the Apollo 17 site. Two additional periods supplied redund.int d.ita fur 
correlation of events occurring at timev of maximum and minimum .ictixity. A 
six-hour period bracketing a solar eclipse was also recorded. Events recorded 

I 
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during the extended listening periods numbered in the thousands, with 1730 events 
recorded that were sufficiently large to warrant detailed studv. Except for 324 
events generated within the lunar module (L.M). and a few possible meteoroid 
impacts, all the events detected appear to be thermal moonquakes. 

The origin of thermal moonquakes has been a mystery since they were first 
observed on the Apollo 1 1 records. At that time, the most likely source was 
thought to be the I.M. Data obtained from the Apollo 14 mission showed that 
there had to be many sources for these events that tend to repeat themselves at 
the same time of each lunation in numbers correlating with the diurnal heating 
cycle (Fig. 1). Because of this repetition they could not be generated by meteoroid 
impacts and because there were multiple sources they could not all be generated 
in the LM, thus the term “thermal moonquakes" (Duennebier, 1973). Because of 
an apparent correlation of activity levels with the topography near each statio;-. it 
was suggested that thermal moonquakes could be generated by down^lope 
movement of the regolith (Duennebier and Sutton. 1974). The calculated source 



SUN ANGLE 

Fig. I. Thermal mooiiquake activity at the Apollo 14. I.S. anj 17 siicv. Ten ile.ereev 
change in sun angle is about 20 hr. .Activity increases abruptly about 2 days after sunrise 
and decreases rapidly after sunset at each station 
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energio imply a change in slope of 1“ per 100 m.y. in a 100 m deep crater (Cooper 
and Kovach, 1975). In this paper it will he shown that steep slopes are not a 
necessary condition for the occurrence of thermal moonquakes but that source 
locations do tend to correlate with crater locations. 

Analysis 

Thermal moonquake data obtained from the Passive Seismic Experiment has 
been discussed by Duennebier and Sutton (1974). A short review of their findings 
is in order. Characteristics of thermal moonquake signals include emergent 
beginnings, durations of one to four minutes, frequencies between 2 and 10 H/, 
and small amplitudes. Signals with nearly identical wave forms repeal each other 
regularly every lunation usually within one hour of the same sun angle. Some 
signals show gradual changes in wave forms from one event to the next and others 
show cyclic changes; some occur once per lunation, others occur many times per 
lunation. About 50 distinctly different types of events were observed at the .Apollo 
14 site (indicating 50 different sources) while more than 250 types are recognized 
at the Apollo 15 site. The short-period component does not operate at the Apt»llo 
12 site, and thermal moonquakes have not been actively studied at the .^pollo 16 
site, although they apparently are visible on the long-period components as well as 
on the short-period component. Four years after the initial study was made at 
Apollo 14, an attempt was made to predict the occurrence of five types of events 
based on when they occurred four years earlier. Three of the five types were 
observed within an hour of when they were predicted. The wave forms had not 
changed noticeably in four years. 

Data from the Apollo 17 LSPE array ha^ been described by Kovach et al. 
(1973) and Cooper and Kovach (1975). These data are an improvement over the 
PSE data in that each event is observed on four seismometers instead of one, thus 
making source location possible. .An excellent set of calibration events at known 
ranges and with known energies are available. Characteristics of thermal moon- 
quakes observed by the LSPE are nearly the same as those observed at the other 
seismic stations. Differences worthy of note are; (1) event activity has a single 
peak near sunset rather than multiple activity peaks earlier in the lunation as at the 
PSE stations (Fig. 1); and (2) while thermal moonquakes at the PSE stations are 
generally restricted to frequencies below 10 Hz, events with energy as high as 
20 Hz are observed at the LSPE (Fig. 2). 


SoilRCt LfK ATION 

The standard technique for seismic source location is to use the relative times 
of phase arrivals at each seismometer and. knowing the velocity structure, 
triangulating to find the source. This method \s not useful fi>r LSPE data because 
there are no clear phases observed; even the largest signals have emergent 
beginnings. Signals observed at one end of the array show no obvious wave form 
correlations with signals from the same event observe^ away, or even 50 m 
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away (see Fig. 2). However, the observed events can be divided into two types, 
those which show constant signal amplitude across the array and those that show 
significant variations in amplitude. Not including events from the LM. about 200 
events fall into the group showing variations, or approximately 15^ of the total 
number of natural events observed. If the rate of change in amplitude with range 
is known, then this amplitude variation can be used (in place of travel time 
differences) to locate events. 

Records from the l-SPE explosive packages (E.P.’s) detonated at ranges from ■ 
100 to 3000 m provide an excellent means of determining the change in amplitude 
with distance. Because high frequencies attenuate faster than lower frequencies, 
care was taken to filter the records sharply (m the frequency domain) at 
frequencies of interest to avoid contamination from energy at other frequencies. 
To eliminate the effects of changes in efficiency of the different E.P.’s in coupling 
into seismic waves, the source energy was not included as a known parameter and 
each E.P. was treated as independent of the others (Nakamura et al.). The 
equation < 

A=krf (I) 



was used as a model w here A is the amplitude at the peak of the signal and r is the 
distance from the source to the seismometer. The value of ct was then found 
which best fit the E.P. data. The resulting analysis shows that the amplitude varies 
inversely as the distance at low frequencies near 5 Hz. (a = - I ). the exponent 
increasing to the square of the distance at 19 Hz (« = —2). The higher frequency 
was chosen for determination of source locations since it is the highest frequency 
at which most nearby events have sufficient energy for measurement and because 
the large variation in amplitude with distance increases the accuracy of source 
locations. The calibration data are shown in Fig. 3. i 

The differences between measured amplitudes for natural events are used to 
^ triangulate using a least squares technique to find the source. Because of the small 
size of the array and a ± I db accuracy in amplitude determination, only those 
i events within about 400 m of the center of the array could be located. As a check 
on the accuracy of the method, 30 events from the LM were also loc ited. About i 
809?^ of these source locations fall within 50 m of the l.M (located IK7 m from the 
center of the array). Thus, while the method hardly gives pinpoint accuracy, it is 
capable of showing regions of high activity. Note that accuracy will decrease as 
distance from the array increases and that the average event size increases as 






Fig. 2. (_;0 Thermal moonquakes recorded by the Apollo 17 I.Sl’F. Reci)rdin«\ of two 
diflerenrhut nearl' <denlical events (177 days apart) are shown a\ recorded by the fi'ur 
LSPF pcophonec Note the Lick of ci'herencc of wave formv acrovs the .trray ar.d the 
high coherence between the two events indicatmtr ttut they b>>th orijjinate at th ; value 
source. Vertical lines are 1 .7 sec apart, (b,) Nearby thermal moonquakev. Note the chanire 
in amplitude and frequency content recorded at each peophone fhese events occurred 
within 100 m of the array, closest to jteophone 4 One event occurred 1 77 days before the 

other. 

/ 
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Fig. 3. Calibration dala for ampliiuJe-distance relationNhip. 7 he data are from the I.SPF: 
evpIoNive package recordings tillered at 2i) H/ The source energy for each event was not 
constrained to allow for changes in etficiency. The least -sviuares curve shown has a slope 

of -2. 


distance increases: i.e.. small events recorded from sources inside the array may 
not be visible when their sources are outside the array. 

The resulting source location map for natural events is shown in Fig. 4. The 
event activity is obviously not random but concentrated in two regu»ns. one inside 
the array and one to the southwest. The concentration inside the irray correlates 
well with the location of “geophone rock,” a 2 m high rock to the southwest of the 
center geophone. However, none of the events from this region are large enough 
to be observed had they occurred more than 100 m from the array, indicating that 
rocks of this size are not the sourc^of strong thermal moonquakes. Comparing 
Fig. 4 with Fig. 5, the concentration of activity to the southeast corresponds 
generally to the location of two subdued craters, appro.ximately 100 and 200 m in 
diameter. The concentration to the northwest cvtrrelates with a crater (Rudolph) 
about 50 m in diameter. Just as important are the features that do not correlate 
with regions of thermal moonquake activity. The areas of astronaut activity in this 
region (concentrated mostly around the LM and the central ,.*.LSEP) station 
slightly to the north of the center geophone) show no signs of any unusual activity, 
thus thermal moonquakes do not appear to be generated by man-made activity. 
The relatively featureless region to the southwest is alstv an area of little activity 
Note again that many events occurring near the edges and i»iitside of the region 
shown in Fig 4 were not located because the change in amplitude acntss the array 
is too small for adequate resolution. Had the array size been I km on a side rather 
than I0() m, many more events could have been located by this method. 
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SoiFRCE Energy 

Previous^eslimatcs of thermal moonquake source energy were hampered by 
lack of adequate calibration events (Duennebier and Sutton, IV74). While this is 
still a problem, the LSFE has the advantage that the calibration events and 
thermal moonquakes are recorded on the same instruments, thus no conversion to 
actual ground motion needs to be made. The calibration events (eight small 
explosions and the Apollo 17 L.M impact) are still not ideal in that the source 
mechanisms for explosion events (including the L.M impact) and for moonquakes 
are vastly different and a much larger fraction of moonquake source energy will 
be converted into seismic waves than for explosives. In addition, no single model 
has yet been found that explains the characteristics v>f the signals. Cooper and 
Kovach (1975) noted that a two dimensional scattering equation failed to predict 
the observed signal rise times beyond a range of about I km. They also show that 
the value of k in Eq. (I) is frequency and range dependent (k decreasing as range 
and frequency increase). While Eq. ( I ) is an adequate model for near-range source 
location, a better model is needed for esl ",ating the source energy of natural 
events. For this purpose the equation 



is used, where A is the amplitude at the peak of the signal envelope. £.j is the 
source energy. / is the frequency. Q is the quality factor, and U is the average 
velocity of the signal peak. 1 he values of /t. £„. r. /. and U are measured and the 
values of k, a. and Q are found by regression. Only data within 1.5 km ((S events) 
are used since beyond this range the value of Lf becomes a complicated function 
of range and frequency. Within 1 .5 km U averages about 25 m/sec. The regression 
to find the values of k. a, and Q was done for two frequency bands (4-13 Hz and 
13-21 Hz) with 120 amplitude measurements used in each (5 frequencies per 
record, 4 records per event, and 6 events). The results yield k = 141 and 86. 
a = -1.02 and -0.97, and Q = 1400 and 740 for the low and high frequency sets 
respectively. 

To obtain source energy estimates for thermal moonquakes. the difference 
between explosives and moonquakes in the fraction of energy transformed into 
seismic waves must be accounted for. It is estimated that thermal moonquakes are 
more efficient at seismic wave generation such that the value of k should be 
increased by a factor of from 32 to 100 (Duennebier and Sutton. 1974). The 
resulting energies, based on the events located earlier and the higher frequency 
band parameters above, range from l()‘ to 5 x 10' ergs per event. Note, however, 
that only the events closest to the array are included, since it is only for these 
events that a reasonably reliable range is available. Events of similar amplitude at 
greater distances are observed implying that thermal moonquakes with source 
energies of 10’ ergs and possibly larger do occur. 

A reliable estimate of the total amount of energy released by thermal 
moonquakes within a given period d(»es not appear possible, since nu>st of the 
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Fib- 4 Source lucaliun map. Density of thermal moonquakc source locations in the region surrouniling the LSPE array arc 
shown. Comparison with Fig. 5 shows a correlation of activity with large craters and with geophonc rock. This map docs not 

include events that are known to originate at the LM. 
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energy is apparently released by events smaller than the detection threshold of the 
system (Duennebier and Sutton. 1974; Cooper and Kovach. 1975). All plots 
showing cumulative number of events vs. (log) amplitude for thermal miH>nquakes 
have slopes steeper than -2. implying that the large numbers of small events 
release more energy than the larger events. As a minimum estimate for total 
energy release from thermal mounquakes. consider the LSPH region where more 
than 100 events of 10’ ergs occur each month within I km\ Extrapolating to the 
whole moon, 10” ergs per year are released, as compared with 10” ergs per year 
for the deep moonquakes (Lammlein el al., 1974). Note that this estimate for 
thermal moonquakes is conservative in that the region around the LSPE site is 
relatively smooth whereas rougher r> ^ions are expected to have more activity. 

For comparison with another energy source, a single 500 gm meteoroid 
impacting at 20 km/sec releases 10” ergs of energy. An impact of this si^e or larger 
is expected about 2000 times per year on the lunar surface (Duennebier et al.. 
1975). This difference in energies implies that the impacting process should be 
dominant in the sculpture of the lunar surface, which is no surprise since the 
surface is covered with craters. However, thermal movement observed as thermal 
moonquakes could be an important process in modification of existing features. 


SouRCK Mechanism 

With the evidence presented in this '^aper and others (Duennebier and Sutton. 
1974; Cooperand Kovach. 1975). the correlation of thermal moonquakes with the 
solar heating cycle appears inescapable. Additional confirmation was obtained 
during a solar eclipse for which the LSPE array was activated. The cooling period 
of the eclipse is characterized by relatively large thermal moonquakes occurring 
at about 25 events per hour (Fig. 6). During reheating after return of sunlight, 
activity is reduced to mostly small events, although activity is still much higher 
than normal. It is difficult to imagine a mechanism other than thermal stress 
variations which could cause such activity. 

As mentioned earlier (and by Duennebier and Sutton. 1974). thermal cracking 
of rocks along zones of weakness can account for only a fraction of the events 
observed. While events are observed from the region around geophone rock, they 
have very small source energies compared to most thermal moonquakes. although 
this is one of the largest rocks in the vicinity. The change in wave form from one 
event to the next during each lunation for some types of thermal moonquakes 
implies movement of the source, thus excluding rocks as sources for these events 
(Duennebier and Sutton. 1974). 

Most thermal moonquakes appear to be generated in the regolith itself 
Duennebier and Sutton (1974) arrived at this conclusion and suggested that the 
dominant source mechanism is triggering of gravitational downslope regolith 
movement by thermal stress variations. While this conclusion still appears ti> be 
true, many of the source locatiims obtained in this p.tper do not correlate with 
steeply sloping surfaces. Apparently the thermal stresses alone supply enough 
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Fi*. 6. Thermal activity at the LSPE during a solar eclipse. Trace amplitude is plotted as 

a function of time in hours. 

energy for generation of thermal moonquakes and gravitational energy is not 
necessary. 

Duennebier (1972) showed that, with reasonable assumptions as to the 
compressibility and strength of the lunar soil, the thermal stresses in the soil 
would be the same order of magnitude as the shear strengths, thus movement of 
the soil can be expected. Where and how these movements occur is a matter for 
debate. One might expect that the lunar soil would expand and contract elastically 
or "fracture” at close enough intervals such that no deformation or movement 
would be noticeable. Yet. if the new values of source energy are correct, using the 
method of Duennebier (1972) the source area will be about 6 cm in radius for a 
thermal moonquake of 10* ergs source energy, and maximum displacement in this 
region will be about 1 mm. This value is reasonable within an order of magnitude if 
one assumes a linear thermal expansion coefficient of .S x 10 *. in w hich case the 
diurnal surface temperature variation of .TSO’K requires 1 mm expansion for every 
50 cm. 

Regions in which one might expect the most movement are where there are 
strong contrasts in temperature or in physical properties; such contrasts exist 
across shadow lines caused by sloping surfaces or other surface features. The 
tendency will be for material on the hot (expanding) side to push into the cold 
(contracting) side. In craters, the shadowed portion is always d«)wnslope of the 
sunlight region such that material will always tend to move downslope (Fig. 7). 
Along sharp crater edges, material from outside the crater will tend to move into 
the crater. Generally, the process will be one of smoothing the lunar surface and 
eradicating discontinuities. Cooper and Kovach (1975) computed a slope degrada- 
tion rate of 1° per 100 m.y. for a 100 m high slope using the method of Duennebier 




Fig. 7. Movement uf the regolith cuu^ed by thermal expanvion acrovs shadovb> The 
motion, shown by arrows, from hot regions to cold (shadowed) regions, will always be 

downslope. 


and Sutton (1974) which assumes that gravitational energy is the only energy 
released. Using this model, fewer events are necessary to change gentle slopes hy 
one degree than steep ones in craters of the same radius since the change in 
gravitational energy per degree becomes less as the slope decreases. Intuitively, 
however, one c.xpects slopes to change rapidly when they are steep and more 
slowly as the erosion continues. It seems likely that gravitational energy release is 
significant only for steeper slopes. For shallower slopes, thermal energy probably 
dominates. If we assume that a typical thermal moonquake displaces 200 cm’ by 
1 mm, then we can calculate how long it will take to degrade a crater by I* by 
computing how much volume is moved in the transition. A 5 m deep crater with a 
28 m radius will have about 130 m’ of material transported an average of 14 m in a 
slope change from 10“ to 9“. This process will require 9x 10'* average thermal 
moonquakes and will take about 9m.y. if 1000 events occur per year. A crater 
with twice the radius (and four times the number of events occurring per year) will 
require four times as long for the same transition. 

CONCt.USIONS 

The movement of the regolith in response to diurnal stresses iy observed as 
small high-frequency seismic events (thermal moonquakes). Most of this move- 
ment is apparently related to regolith materials moving downslope as sunlit areas 
expand and shadowed areas contract. Sources of thermal moonquakes located 



near the Apollo 17 LSPE array suggest that slopes at angles much less than the 
angle of repose can be sources for these events indicating that thermal energy, 
rather than gravitational energy, is the main driving force. Ihe effect of this 
activity over long periods of time ill be to fill craters, lessen slopes, and generally 
eradicate sharp discontinuities. 

The importance of this mechanism relative to others suggested for lunar 
erosion is hard to determine. Other mechanisms include impact cratering (Arnold. 
I97.S). seismic shaking (Houston et ai, IV72; Schultz and Gault. I97.S). and 
electrostatic transport (Gold. 1%0: Criswell. 1972). Downslope movement is well 
documented by photographs and astronaut observations (Mattingly. 1973). and the 
fact that fines move downslope faster than large rocks can be seen from the 
observation that fillets occur preferentially on the uphill side of rocks 
(Muehlberger et at., 1973). Crater life times are also significantly shortened on 
slopes (Basilevsky. 1976). In an earlier paper (Duennebier and Sutton. 19M) it was 
suggested that more than enougn energy was available in the thermal moonquake 
mechanism to account for significant erosion of lunar slopes. .Although the source 
energy estimates have been lowered significantly since that paper, this conclusion 
still appears to be valid. 
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